Velocity renormalization of nodal quasiparticles in <i-wave superconductors 
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Gapless nodal quasiparticles emerge at a low-energy regime of high-T c cuprate superconductors 
due to the d x 2_ y 2 gap symmetry. We study the unusual renormalizations of the Fermi velocity vf and 
gap velocity da of these quasiparticles close to various quantum critical points in a superconducting 
dome. Special attention is paid to the behavior of the velocity ratio, va/vf, since it determines 
a number of observable quantities. We perform a renormalization-group analysis and show that 
the velocity ratio may vanish, approach unity, or diverge at different quantum critical points. The 
corresponding superfluid densities and critical temperatures are suppressed, slightly increased, or 
significantly enhanced. The effects of three types of static disorders, namely, random mass, random 
gauge potential, and random chemical potential, on the stability of the system are also addressed. 
An analogous analysis reveals that both random mass and random gauge potential are irrelevant. 
This implies that these fixed points of the velocity ratio are stable, and hence observable effects 
ignited by them are unchanged. However, the random chemical potential is marginal. As a result, 
these fixed points are broken, and thus, the instabilities of quantum phase transitions are triggered. 



PACS numbers: 73.43.Nq, 74.72.-h, 74.25.Dw 



I. INTRODUCTION 

It has been well-established that high-T c cuprate su- 
perconductors have an anisotropic d x 2_ y 2 energy gap. 
This gap vanishes at four nodes (±|,±~), so the low- 
energy elementary excitations are gapless nodal quasi- 
particles (QPs) with a linear energy spectrum. These 
nodal QPs are responsible for many low-temperature 
thermal and transport properties of the superconducting 
phase of high-T c superconductors For instance, the 
specific heat exhibits a linear temperature dependence, 
C(T) oc T 2], which is quite different from that of con- 
ventional s-wave superconductors, and indeed has played 
a crucial role in the determination of d x 2_ y 2 gap symme- 
try. In the superconducting state, nodal QPs are revealed 
by numerous experiments, especially angle-resolved pho- 
toemission spectroscopy (ARPES) Q and transport mea- 
surements 0, [H, to be well-defined, with extraordinary 
long mean- free paths [l[ . A residual short-range interac- 
tion between nodal QPs leads only to a quite weak damp- 
ing rate, oc max (w 3 , T 3 ) , and insignificant corrections to 
fermion velocities [6|. 

In 1999, ARPES measurements by Valla et al. re- 
vealed a fermion damping rate, oc T, in the optimally 
doped cuprate superconductor P^Si^CaC^Os+a, which 
indicates a marginal Fermi liquid behavior [8] . This un- 
expected finding is apparently in sharp contradiction to 
theoretical predictions. Such a strong damping can only 
be caused by some kinds of soft (gapless) bosonic modes 
Q. In order to explain this experimental finding, Vo- 
jta et al. proposed that the soft boson may arise 
from a quantum phase transition from a pure d x 2_ y 2 su- 
perconducting state to a new X superconducting state. 
Generically, there are a number of candidates for the pro- 
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posed X state, and it seems difficult to uniquely deter- 
mine which is the correct one. Vojta et al. }9Hll| con- 
sidered all possibilities for the order parameter of the X 
state, and ruled out most of them after carrying out care- 
ful symmetry analysis and field-theoretic calculations. 

Nodal QPs have two velocities: the Fermi velocity vf 
and the gap velocity va [Hj]. Extensive experiments 
have determined that they are not equal to each other, 
VFjt va- Indeed, transport and ARPES measurements 
[JlH found that the velocity ratio, va/vf ~ 0.1, in most 
high-r c superconductors. Remarkably, a number of im- 
portant observable quantities depend on such a velocity 
ratio pi, including the superfluid density [l3j], critical 
temperature T c [13[ , and electric and thermal conductiv- 
ities [HI, [lj, [lfij]. Any unusual renormalization of this 
velocity ratio will give rise to considerable changes in 
these observable quantities. In the presence of transition 
from a pure d-wave superconducting state to a new X su- 
perconducting state, the fluctuation of the new order pa- 
rameter can forcefully couple to gapless nodal QPs near 
the quantum critical point, which may lead to nontrivial 
velocity renormalizations. If this ratio deviates strongly 
from its bare value, these physical quantities will be ei- 
ther enhanced or suppressed. 

In this paper, we focus on the unusual renormalizations 
of fermion velocities caused by the critical fluctuations of 
different X order parameters in the d-wave superconduc- 
tor. We are particularly interested in the low-energy be- 
havior of the velocity ratio va/vf- Inspired by the study 
of velocity renormalization in graphene [16l4l8| . it can 
be analyzed by the renormalization-group (RG) method 
p^ - [2l| . The unusual ratio va/vf will have significant 
impacts on transport properties. Moreover, the trans- 
port properties of nodal QPs are largely determined by 
scattering due to impurities. It is, therefore, necessary 
to examine the influence of various impurity potentials, 
other than the interaction between nodal QPs and X or- 
dering, on the behavior of fermion velocities. Based on 
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the coupling between nodal QPs and disorders, there are 
three types of disorders in a <i-wave superconductor: ran- 
dom masSjjandom gauge potential, and random chemical 
potential [22|]. We include these disordered potentials in 
our model and explore their effects by means of the RG 
method. 

Clean-limit systems are considered first. We employ 
three types of vertex matrices, M = t x , t v , and r z , 
to denote various quantum phase transitions. [Quantum 
phase transitions from a pure d x 2_ y 2 to a new X super- 
conducting state can be classified by vertex matrices, M, 
between nodal QPs and X order parameters [HI], shown 
in Eq. and henceforth we dub the new X super- 

conducting state the X matrix state.] By implementing 
the RG analysis [l9l - [2"lj |. we obtain three distinct fixed 
points of the velocity ratio due to the fluctuations be- 
tween nodal QPs and X order parameters in the vicinity 
of the quantum critical points, which are va/vf — > 
20], va/vf — > 1, and va/vf — > °o for M = t x , t v , and 
r z , respectively. The case va/vf — > has been discussed 
in recent publications [2(J [H| , the others will be focused 
on in this paper. Since many physical properties display 
va/vf dependence, observable effects kindled by these 
interesting fixed points should be expected. Indeed, we 
find that superfluid density and critical temperature are 
sensitively influenced by approaching these fixed points. 
Both are suppressed, slightly increased and significantly 
enhanced for cases M — t x , t v , and r z , respectively. 

Whether the results in the clean limit are stable against 
the disorders are also investigated. In practice, the disor- 
der effects induced by various kinds of scattering are in- 
evitably present in the low-temperature transport prop- 
erties of an interacting electron system. If the fixed 
points are changed or even broken by the disorder effects, 
the behavior of physical observables will be considerably 
affected. Hence, it is imperative to examine the disorder 
effects on RG flows of fcrmion velocities. In general, there 
are three types of disorders (random mass, random gauge 
potential, and random chemical potential [22j] ) coupled 
to gapless nodal QPs in the d-wave superconductor. The 
impacts of these disorders on the low-temperature trans- 
port properties of nodal QPs have been studied exten- 
sively [23, [24| . In our case, the RG flows of fermion veloc- 
ities will be influenced by these disordered potentials, and 
in the meanwhile, the RG flows of strength parameters 
of fermion-disorder couplings are determined by fermion 
velocities. Therefore, we should self-consistently com- 
pute the flows of fermion velocities with disorder strength 
parameters. After a detailed RG analysis of the inter- 
play between X order parameter fluctuations and disor- 
der scattering, a series of coupled RG equations of Fermi 
velocity vp , gap velocity wa, and disorder strength pa- 
rameter vr are derived. Based on numerical calculations, 
we learn that both random mass and random gauge po- 
tential are irrelevant. This signifies that these two types 
of disorders do not change the flows of fermion velocities. 
Accordingly, the corresponding fixed points and physical 
observables are stable. However, the random chemical 



potential is marginal. Consequently, these fixed points 
will be destroyed and hence the instabilities of quantum 
phase transitions are signaled. 

The rest of the paper is organized as follows. The ef- 
fective field theory and the corresponding Feynman rules 
are presented in Sec. |H] We calculate the self-energy and 
vertex corrections in Sec. IIIII and Sec. IIV1 respectively. 
A detailed RG analysis is given in Sec. E] which is fol- 
lowed by discussions of the numerical solutions of the RG 
equations in Sec. I VII and of the behaviors of the super- 
fluid density and critical temperature caused by velocity 
renormalization in Sec. IVHI Finally, we briefly summa- 
rize our results in Sec. IVHT1 



II. EFFECTIVE FIELD THEORY OF 
QUANTUM CRITICAL PHENOMENA 



We begin with the action 



S — S'vp + S,h + S< 
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where the free action for nodal QPs is 
Sv = [ T^^Tr^^iai-iu + v F k x r z +V A kyT X )^i a 



(2tt) 2 2tt 
(2tt) 2 2^ 



^l a {-iu + v F k y T z + v A k x T x )^ 2a , (2) 



where r^'^ 2 ) denote Pauli matrices. The linear disper- 
sion of Dirac fermions originates from the d x 2_ y 2-w&ve 
symmetry of the energy gap of cuprate superconductor. 
Here, the Nambu spinor ^\ represents nodal QPs ex- 
cited from the (§,§) and (—§,—§) nodal points, and 

^>2 represents the other two nodal points; w is a Matsub- 
ara frequency in the zero-temperature limit, k XiV describe 
the wave vector from the nodal points and have been ro- 
tated by \ and Vf and va are the Fermi velocity and 
the gap velocity, respectively [20[ • The repeated spin in- 
dex a is summed from 1 to Nt, the number of fermion 
spin components. The ratio va/vf ~ 1/20 between the 
Fermi velocity and the gap velocity is determined by ex- 
periments pun. 

With the help of group-theoretic classification, Vojta 
et al. [9HT1I] pointed out seven possible quantum phase 
transitions from a pure d x 2_ y 2 superconducting state to 
X matrix states. The effective action S$, which describes 
the X order parameter in real space is 



S,i 



d 2 xdr 



-(V0) 2 



24 



4 },(3) 



where t is imaginary time and c is the velocity. The mass 
parameter r drives the system to undergo quantum phase 
transitions, with r = defining the zero-temperature 
quantum critical point, uq is the quartic self-interaction 
strength. The interaction between nodal QPs and 
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FIG. 1: The polarization function for the X order parameter. 
The solid line represents the fermion propagator and the wavy 
line represents the boson propagator. 




(a) (b) 

FIG. 2: One-loop fermion self-energy correction due to (a) 
X-order-parameter fluctuation and (b) disorder. The dashed 
line represents disorder scattering. 



the X order parameter <j> is described by a Yukawa cou- 

piing (mm 

5^ = / d 2 xdr |a (K*i^i*i + ^lM 2 ^ 2 )} , (4) 
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with Ao the coupling constant. There are seven possible 
X matrix states, which can be distinguished by different 
choices of matrices Mi and M 2 [l(| E| ■ It was shown in 
Refs. Ql [Hi that two states out of these seven candidates 
are irrelevant. Therefore, we only need to concentrate on 



Mo 



= T y- 

»; (iv) 



the remainder of the states: (i) M\ 
(ii) M 1 = t v , M 2 = -t v ; (iii)Mi = t x , M 2 
Mr = t z , M 2 = -t z ; and (v)Mi = t x , M 2 = -r x . 

The state with Mi = M 2 = t x corresponds to a ne- 
matic state, which has been extensively investigated in 
a number of papers [2(| HH [UHll] . The critical fluctua- 
tion of such nematic order leads to an extreme velocity 
anisotropy and other unusual properties. The influence 
of disorders was studied in Ref. 21]. In this paper, we 
consider the other four states. Fortunately, when calcu- 
lating the fermion self-energy and boson polarization, M\ 
or M 2 always appears in pairs. Hence the contributions 
do not depend on the signs of Pauli matrices. For in- 
stance, cases i and ii or iii and v would lead to the same 
results. In order to simplify the discussion, we can ignore 
the possible minus of M 2 to obtain a compact classifica- 
tion. Therefore, we have only two possibilities: M = t v 
and M = t z . (The basic conclusions are independent of 
this simplification.) 

To analyze the above field theory, we consider a general 
and large fermion flavor Nf, and apply a 1/Nf expansion. 
The free fermion propagator is 



G°(k, w ) 



1 



— iuj + vpk x r z + v/±k v T x 



(5) 



for nodal QPs (the free propagator for nodal QPs ^ 2 
can be written similarly). 

To the leading order of 1/Nf expansion, the polariza- 
tion function is shown in Fig. ([T]) and symbolizes the 
integral 



n(q,e) 



d 2 k duj 



Tr[MG% (k, w)MG% (k + q, u + e)] , 



(2tt) 2 2tt 

After straightforward calculations [2(| Hl| , we have 
n(q,e) = n»(q,e) 



16vfva 
I 



€ z + v 2 F ql + 



1&vfv\ 

for the case of M = t v , and 



(6) 



n(q,e) 



n z (q, e ) 
1 



I&vfva (e 2 + vpql + v\q^) x l 2 
1 (c 2 - 



vWx) 



l&v F VA (e 2 + v F q% + v\q£) x l 2 



(7) 



for M = t z . In a low energy regime, the polarization 
function is linear in \q\ and, therefore dominates over 
the g 2 -term. Near the quantum critical point, we keep 
only the mass term and assume that (f> — > 4>/Xo and 
r — ► Nfr\ 2 , leading to 

S = Sv+Jd 2 xd T {^(f> 2 + <f>[i&l a M*i a + *t a MV 2a }}.(8) 

After integrating out the fermion degrees of freedom, the 
effective action for the scalar field (X order parameter) 
becomes 



Nf~ 



d :i 



(J 



(2nf 



[r + n(q)M q )\ 2 +0( 



(9) 



Now the effective propagator of the X order parameter 
is 



G<j>(q, e) 



1 



n(q,e) 



(10) 



at the quantum critical point r = 0. Propagators (JSJ) and 
(jTUJ) will be utilized in the following calculations of the 
fermion self-energy and RG equations. 

In almost all realistic condensed matter systems, disor- 
ders are present and play significant roles in determining 
the low-temperature behaviors. In the present problem, 
the nodal QPs can interact with three types of random 
potentials, which represent different disorder scattering 
processes. According to the coupling between nodal QPs 
and disorders, there are three kinds of random potentials 
in d-wave superconductors: random mass, random chem- 
ical potential, and random gauge potential. All these 
types of disorders have been investigated in the contexts 
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of the d-w&ve cuprate superconductor [23|, |24j , quantum 
Hall effect , and graphene [H, [5(| ■ In the subsequent 
analysis, we consider all three sorts of disorders. 

The coupling term between the fermion field and a 
random field A(x) can be written as [22[ 



/ 



ci 2 x* t (x)rV'(x)A(x) 



(11) 



The matrixes T are r = I, T — t v , and T = (t x : t z ) for 
a random chemical potential, a random mass, and a ran- 
dom gauge potential, respectively. We assume that the 
random potential A(x) is a quenched, Gaussian white- 
noise field with the correlation functions 

(A(x)> = 0; (A(x 1 )A(x 2 )) = gv 2 8 2 ( Xl - x 2 ), (12) 

where the dimensionless parameter g represents the con- 
centration of impurity, and the parameter vr measures 
the strength of a single impurity. It will be convenient 
to redefine the random potential as A(x) — > vrA(x), and 
then write the fermion-disorder interaction term as [22j 



S dis =v r J d 2 x*t( x )rV,(x)A(x), 



(13) 



with the random potential distribution 

(A(x)) = 0; (A(xi)A(x 2 )> = gS 2 (yn - x 2 ). (14) 

Then by studying the vertex correction to the fermion- 
disorder interaction term, we can obtain the RG flow of 
disorder strength. After a Fourier transformation, the 
corresponding action for the fermion-disorder interaction 
has the form 

Sdis = vtJ d 2 kd 2 kidw* t (k,w)r*(ki,w)A(k- ki).(15) 



This action is analyzed together with the actions ([2]) , ([8j , 
and (j^J). In order to carry out perturbative expansion, 
both g and Vr are assumed to be small in magnitude, 
corresponding to the weak disorder case. 

Before making RG analysis of the velocities and dis- 
order strength parameter, we calculate the one-loop 
fermion self-energy and vertex corrections in the follow- 
ing two sections. 



III. FERMION SELF-ENERGY CORRECTIONS 

The interplay between X-order-parameter fluctuation 
and random potentials can yield self-energy corrections 
to the free propagator of nodal QPs, which are described 
by the Dyson equation 

(k,u) = -iuj + v F k x r z + VAk y T x 

-Ex(k,w)-E d i.(k,w), (16) 

where self-energy functions £x(k, w) and £dis(k, uj) come 
from X order parameter fluctuation and disorder scatter- 
ing, respectively. To the leading order, the corresponding 
Feynman diagrams of self-energy are presented in Fig. 
©• 

Employing the method of Ref. [13] , we obtain 



cE x (k,w) 



dlnA 



Ci(-iw) + C 2 v F k x r z + C 3 v A kyT x , (17) 



where 



for M = r» and 



where 



n 2(v A /v F ) /"» [ 2 \ a -x 2 + cos 2 9 - (v A /v F ) 2 sin 2 9 

C* 2 = 3 / dx I dO—— — — — ' , 2oX2 G(x,6), (19) 

N f n J-oo Jo {x 2 +cos 2 9+(v A /v F ) 2 sm z 9) 2 



Nfir 3 J ( x 2 + cos 2 9 + (va/v f ) 2 sin 2 9) 2 V h V ' 

\jx 2 + cos 2 9 + (v A /v F ) 2 sin 2 + \jx 2 + sin 2 9 + (v A /v F ) 2 cos 2 0, (21) 



d ?i n k A^ ) = Dl (~ iuJ ) + D2V F k * rZ + D 3 v A k y r x , (22) 



Dl = 2 ^y) r dx r d ^-cos 2 ^-(W^) 2 sin^ 

N f n* J (x 2 +cos 2 9+(v A /v F ) 2 sm 2 0) 2 ^ h V ' 
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for M 



D 2 
D 3 



2(v F /v A ) 

NfTT 3 

2{v F /v A ) 

NfTT 3 



,-2n 



dx 



de 



— oo 
oo 
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x 2 - cos 2 9 + (va/vf) 2 sin 2 9 
(x 2 + cos 2 9 + (v A /v F ) 2 sin 2 6») : 



■J*(x, 



dx / <20 

o JO 

I- sin 2 9 



— X i-ua v -j- \u A / up j fin u 

(x 2 + cos 2 + {v A /v F ) 2 sin 2 0) 2 (X ' j ' 
x 2 + cos 2 6> 

2 cos 2 



cos 2 + (v A /v F ) 2 sin 2 



^x 2 + cos 2 9 + {v A /v F ) 2 sin 2 9 \j x 2 + sin 2 + (v A /v F )'< 



(24) 
(25) 
(26) 



The fermion self-energy due to disorder £di S («w) can 
be computed as 



£ dis (zw) = gv 2 



J (0 rG °( k ' w ) r 



2irv F v A 



i In A. 



(27) 



According to this result, we infer that £dis(*w) exhibits 
the same behavior for all possible choices of T and is ac- 
tually Af -matrix independent. Another striking feature 
is that Sdis(jw) does not depend on momentum, which 
reflects the fact that quenched disorders are static. This 
leads to 



dT, dis (iu>) . 



where 



C a 



2ttv f v a 



(28) 



(29) 



IV. VERTEX CORRECTIONS 

The fermion-disorder interaction parameter vr is also 
subjected to RG flow. To obtain its flow equation, we 
need to calculate the fermion-disorder vertex corrections. 
Formally, the vertex correction has the form 



v r T' = v r T + T i + V i , 



(30) 



where I\ represents the vertex correction due to X or- 
der parameter fluctuation and Vi represents the vertex 
correction due to disorder interaction, and the quantities 
denoted i = 1,2,3 correspond to the random chemical 
potential, random mass and random gauge potential, re- 
spectively. The Feynman diagrams are shown in Fig. 
([3]). Both cases, M — t v and M — t z , are calculated 
explicitly in the following for all three types of disorders. 



A. Random chemical potential 

To compute the vertex correction owing to X ordering, 
we take advantage of the method proposed by Huh and 



Sachdev [20j. At zero external momenta and frequencies, 
the vertex correction is expressed as 



Fr = vr I ^H{Q)^ 1 q 



2 ' 

A 2 ) ' 



(31) 



where JC(x) is an arbitrary function with /C(0) = 1, and 
it falls off rapidly with x, e.g., JC(x) — e~ x . However, the 
results are independent of the particular choices of JC(x). 
There is a useful formula (20| . 

where 



H(Q) = ~A1~ T 



Nf (— iv F x + v F cos9t z + v A sin St 1 



(— ivpx + vp cos9t z + v A siti9t x ) Tl(Q) 

(33) 

where i = y and z denote types M = t v and M = t z , 
respectively. Here, matrix I corresponds to the coupling 
between nodal QPs and the random chemical potential. 
It will be replaced by t v in the case of random mass 
and t x,z in the case of random gauge potential. After 
straightforward calculation, we have 



where 
C 5 



C 5 v r l, M = r y 
D 5 v r l, M — t z . 



/•27T 

dx I d9 

v J 



and 



D* 



dT x 
dlnA 



2(v a /vf) 

NfTT 3 

(x 2 — cos 2 9 — (v a /vf) 2 sin z 
% (x 2 + cos 2 9 + (v A /v F ) 2 sin 2 9) 2 
-Ci, 

2 ^ r dx r d9 



(34) 



G(x,9) 



(35) 



NfTT 3 Jo 



( [x 2 - cos 2 9 - (v A /v F ) 2 sin 2 9) 
* (x 2 + cos 2 9 + (v A /v F ) 2 sin 2 9) 2 
-D x . 



T(x,9) 



(36) 
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(a) 



(b) 



FIG. 3: Fermion-disorder vertex correction due to (a) the X 
order parameter and (b) the disorder parameter. 



The vertex correction due to averaging over disorder is 

V^gv 2 J ^IG%(cj,p)v T IGl(uj,p + k)I. (37) 

Again, one should substitute a certain Pauli matrix for 
matrix I in the case of random mass or random gauge 
potential. Taking the external momentum k = and 
keeping only the leading divergent term, we have 



dV 1 
rflnA 



where 



for both types. 



C r = v ll = c 



(38) 



(39) 



B. Random mass 

Calculation of the vertex correction in the case of ran- 
dom mass parallels the process presented above, so we 
just state the final result. The A-ordcring-induced ver- 
tex correction is 



dY 2 
rflnA 



C 6 v r T v , M = t v 
D 6 vrTy, M = t z . 



(40) 



where 
C 6 



do p2tt 

dx / d6 

co JO 

2 a - (va/vf) 2 sin 2 9) 

* (x 2 + cos 2 9 + (v A /v F ) 2 sin 2 9) 2 
C3 — C\ — C2, 



2(v A /v F ) 

A/7T 3 

(-X 2 



COS 



G(x,0) 



(41) 



The disorder-induced vertex correction, in both cases, 
can be written as 



dV 2 
dlnA 



where 



vh = c 

2-KVpVA 



(43) 



(44) 



C. Random gauge potential 

The random gauge potential has two components, 
characterized by t x and t z respectively. For the t x com- 
ponent, the A-ordering contribution to vertex correction 
is 



dT 3 
d\nA 




(45) 



where 



/CO />27T 
dx I d6 
-co Jo 



2(v A /v F ) 
Nfir 3 

(x 2 + cos 2 9 - (va/vf) 2 sin 2 1 
(x 2 + cos 2 9 + (va/vf) 2 sin 2 9) 2 ' 
-C 3 , (46) 



D 4A = 



2(vf/va) 



/co r-2 
dx \ 
-co Jo 



d9 



Nfir 3 

(x 2 + cos 2 9 - (va/vf) 2 sin 2 9) 
(x 2 + cos 2 9 + (va/vf) 2 sin 2 9) 2 
= -D 3 . 

For the t z component, we have 

dT 1 C iB vvT z , M 



T(x,9) 



rflnA 



D 4B vtt z 



M = t z . 



(47) 



(48) 



where 
C4B 



r 2TT 



d9 



(x 2 - cos 2 9 + (va/vf) 2 sin 2 9) 



G(x,9) 



(x 2 + cos 2 9 + (va/vf) 2 sin 2 9) 2 ' 
-C 2 , (49) 



D 6 = 



1" J-oo Jo 



d9 



2(v f / va) 
N f 

(x 2 + cos 2 9 + (va/vf) 2 sin 2 9) 
X (x 2 + cos 2 9 + (v A /v F ) 2 sm 2 9) 2 (X ' ' 
D 2 -D 1 -D ?> . 



(42) 



2(v f / va) 
N f 

(-X 2 + cos 2 



1^1 dx de 

J-00 Jo 

(va/vf) 2 sin 2 t 



F(x, 



(x 2 + cos 2 9 + (va/vf) 2 sin 2 9) 2 
= -D 2 . (50) 
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The disorder contribution can be calculated similarly. 
For both the t x and the t z components, we have 



V 3 (oj) — finite, 
for types M = t v and M = r z . So 
dV 3 



din A 



0. 



(51) 



(52) 



RG ANALYSIS 



In this section, we make an RG analysis of the fermion 
velocities and disorder strength and then derive RG equa- 
tions. To this end, it is convenient to perform the scaling 
transformations [19l - [2l| 



k 



k'e 



V 1<2 (k,oj) 
0(k,w) 



*i i2 (k / ,« / )e» / o ,( * _, " ),H , 
0'(k',o;')e' J o ! (5-w)<« j 



(53) 
(54) 

(55) 
(56) 



where 6 = e~ l with Z > 0. The parameters r)f and rji, are 
determined by the self-energy and X-ordering-fermion 
vertex corrections. Note that the energy is required to 
rescale in the same way as the momentum, so the fermion 
velocities are forced to flow under RG transformations. 

The standard procedure for assigning the scaling trans- 
formation of a field operator when the energy and mo- 
mentum are rescaled, according to the spirit of RG theory 
19], is to keep its kinetic term invariant. Since the ran- 
dom potential A(x) does not possess its own kinetic term, 
it actually does not work in the present problem. In or- 
der to find out its scaling behavior, we write the Gaussian 
white-noise distribution in the momentum space as 



(A(k 1 )A(k 2 )) = ff «5 2 (k 1 +k 2 ). 



(57) 



When the momentum k becomes &k, the delta function 
is rescaled to 

<5 2 (kx + k 2 ) ->■ <5 2 (6k! + bk 2 ) = &- 2 £ 2 (k! + k 2 ). (58) 



If we require that the disorder distribution, Eq. (|57p . is 
invariant under scaling transformations, then the random 
potential should be transformed as 

A(k) -> & _1 A(k). (59) 

Now we have to assume that 

A(k) - A'(k')e l . (60) 

In the light of the RG technique introduced in Refs. 
fl6| . [Isj . and [l||, the momentum shell between bA and 
A will be integrated out, while keeping the — iu) term in- 
variant. From the message of type M — t v ordering and 



disorder contributions to the fermion self-energy func- 
tion, we have 



bA 



d 2 kdu& 

r bA 



—iio — Ci(—ioj) In — + C a (—iio) In — 
6A s 



bA 



/ d 2 kdu;¥ + {C g - Ci)Z]* 

/bA 
d 2 kdLJ^(-ioj)e( C *- Cl '> l t>. 



(61) 



After the scaling transformation, this term should go 
back to the free form, so that 



vf = c g - a 



(62) 



The kinetic terms should also remain invariant under 
scaling transformations, which leads to 



dvp 
~df 
dvA 
~df 



(C\~C 2 -C g )v Fl 
(Ci-C 3 -C> A . 



(63) 
(64) 



Based on these expressions, the ratio between the gap 
velocity and the Fermi velocity is given by 



d(v A /v F ) 
dl 



(C 2 -C 3 )(v A /v F ). 



(65) 



By replacing C with D in Eqs. (|5T| [6"2l 1551 R%1 [55) . we 
could get similar equations for the case M = t z . 

The disorder strength parameter vr enters the above 
expressions. Because of the interplay of X ordering and 
disorder, this parameter also runs under RG transforma- 
tions. The flow equation depends on the type of disorder, 
which is studied in the following. 

We first consider the case of the random chemical po- 
tential for M = t v . The bare fcrmion-disorder action 



J d 2 kd 2 k 1 dw* t (k,a;)r*(k 1 ,w)A(k-k 1 ). (66) 



Taking into account corrections due to ordering and dis- 
order interactions yields 



bA 



(^k^kidwtf^k, 



vrl — Cifrlln — 
bA 



CgVrlhi-^- 
a bA 

r bA 



*(ki,cj)A(k-ki) 

/Oil 
d 2 kd 2 kiduj&(k,cj)vrl[l + {C g - d)l] 

x*(ki,w)A(k- ki) 

rbA 

I d 2 kd 2 kidw*t(k,a;)urIe (C9 " Cl)i 
x*(ki,w)A(k-ki). (67) 
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where 




FIG. 4: (Color online) Flows of vf,va, and va/vf at a rep- 
resentative initial value vao/vfo = 0.1 for M = t v . The 
conclusions are independent of these concrete values. 



ZTTVpVA 



i = 1,2. 



(75) 



Here, the equations denoted by i = 1,2 correspond to 
the t x and t z components of the random gauge poten- 
tial, respectively. Their M — r z -case counterparts are 
conveniently obtained by substituting D for C in Eqs. 

ClD [H [73 M E3- 



VI. NUMERICAL RESULTS 

In this section, the numerical solutions of the RG equa- 
tions are presented and the physical implications of these 
results are also discussed. We first consider the clean 
limit and then include random potentials. 



After redefining the energy, momentum, and field opera- 
tors, we are left with 

,A 

/ d 2 \s!<P\s!idu)'^{\s!,w')v v l 

xe^ c "- Cl ^'(\s' 1 ,Lj')e-^ l A'(\s' - k[). (68) 

Since 77/ = C g — Ci, it is easy to obtain the following RG 
flow equation for v^, 



= 



(69) 



Evidently, the parameter vr does not flow and thus can 
be simply taken to be a constant. By applying the similar 
steps, we could get the same result for M = t z . 

In the case of random mass, the running equations 
for fermion velocities have the same forms as Eq. (|63[) 
and Eq. (I64[) . However, the flow equation for disorder 
strength parameter is different f rom Eq. (1551) , and would 
be recast into 



dvr 
~df 



(C 3 - C 2 - 2C g )vr, 



M 



(D 2 -D 3 -2D g )v r , M = t* 



(70) 



which couples self-consistently to the flow equations of 
the fermion velocities. 

Following the steps presented above, we can derive the 
corresponding RG equations in the case of random gauge 
potential for M = t v 



dvp 
~df 
dvA 
~df 



= {Cx-Cz-C^vf, 

= (Cx-Cs-C^VA, 



(71) 
(72) 



which couple to the flow equations of disorder strength 
dvri 



dl 
dv T 2 
~df 



[{d - C gl ) - Cyvn, (73) 
[(Ci - C g2 ) - C 2 ]vr2, (74) 



A. Clean limit 

By analyzing the coupled RG equations of fermion ve- 
locities vf and ua introduced in Sec. IVl with clean limit 
g = 0, we can obtain the running flows with decreas- 
ing energy scale, i.e., growing scale of /. An extreme 
anisotropy of fermion velocities, (va/vf)* — 0, caused 
by the nematic order parameter (M = t x ), was found in 
Ref. [20( 1 . We subsequently list the results of two other 
cases. 

In the case of M = r v , the velocity ratio va/vf flows to 
a new fixed point, (va/vf)* = 1, at the lowest energy, as 
shown in Fig. This implies that the system becomes 
isotropic in the low-energy regime. 

Regarding M — t z , we find another extreme fixed 
point with vf/va in the low-energy regime, as de- 
picted in Fig. ([5]). (The ratio Vf/va decreases rather 
slowly as I increases, and this is discussed in lVI Bl ) This 
indicates that the A-ordering quantum phase transition 
in the d-wave superconductors is accompanied by the ap- 
pearance of an infinite velocity anisotropy. It is inter- 
esting to compare this extreme anisotropy with that of 
va/vf — > driven by the critical nematic fluctuation 



B. Including disorder 

The conclusions in the previous subsection are valid for 
clean systems. In fact, disorders are present in almost all 
realistic condensed matter systems and play important 
roles in determining the low-temperature behaviors. In 
the current problem, the nodal QPs can interact with 
three types of disordered potentials as presented in Sec. 
HT1 which represent different disorder scattering processes. 
In the general analysis that follows, we consider the in- 
fluence of all these kinds of disorders. 

By paralleling the clean limit, the coupled RG equa- 
tions consisting of the fermion velocities vf, and ^a, and 
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FIG. 5: (Color online) Flows of Vf,va, and va/vf at the 
representative initial value vao/vfo = 0.1 for M = t z . The 
conclusions are independent of these concrete values. Inset: 
Flow of vf/va for large / (> 600). 
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FIG. 6: (Color online) Flows of Vf,va,va/vf, and vr in the 
presence of random mass at two representative initial values, 
vao/vfo = 0.1, and g = 10 -3 , for M = t v . The conclusions 
are independent of these concrete values. Inset: Flow of vr 
for large I (> 600). 



disorder strength parameter vr derived in Sec. |V] can 
be numerically solved. The disorder effects at a nematic 
quantum critical point (M — t x ) has been studied in Ref. 
|21j . In the following, the corresponding effects for the 
cases M = t v and M = t z are addressed. In the present 
problem, the combining factor gv^, represents the disor- 
der strength, not only vr, and our analysis is valid for 
weak disorders with small gv"^, and thus both g and vr 
are assumed to be small in magnitude. To obtain the 
compact plots, we have measured vr with vro in Figs. 
©, © and ©. 

The random mass is considered first. The numeri- 
cal calculation is proven to be costly and the decreas- 
ing rate of vr is small in the case of M — t v as pre- 
sented in Fig. (j6]). Although the flows are shown in finite 
I, the tendency is straightforward. The approximately 
analytical discussion also supports this, which includes 
vr(l)/vr « 0.191132 exp{-3.65314 x 1CT 5 (113726 + I)} 
(I > 600) at a representative value g = 10~ 3 . There- 
fore, vr flows to when I approaches infinity. As can be 
easily seen from Fig. J7]) for the case M — t z , vr de- 
creases more rapidly. Despite flowing a little slower, the 
fermion velocity ratio vf/va eventually vanishes in the 
low energy limit. 

We next discuss the case of the random gauge poten- 
tial. By carrying out analogous steps, we come to the 
similar conclusion that the random gauge potential can- 
not qualitatively change the running behavior of Vf, and 
va for both M = t v and M = r z , which are depicted in 
Fig. ©. 

Learning from Eq. (|69p . we are informed that the dis- 
order strength parameter vr, unlike the random mass 
and random gauge potential, is marginal in the presence 
of a random chemical potential. (Since this property is 
independent of the types of matrix states, the indices de- 
noting types of matrix states are discarded without loss 
of generality in the following). This indicates that vr 
does not flow as I grows and hence should be kept as a 



1.5 




FIG. 7: (Color online) Flows of vf,va,va/vf, and vr in the 
presence of random mass at two representative initial values, 
vao/vfo = 0.1, and g — 10~ 3 , for M = r z . The conclu- 
sions are independent of these concrete values. Inset: Flow of 
vf/va for large I (> 600). 



constant. Consequently, the influence of scattering due 
to the random chemical potential can not be neglected. 
As expected, the flows of velocities vf and va are heav- 
ily dependent on the magnitude of vr- According to ex- 
pression (|63|) . it seems that the running behavior of the 
velocity ratio va/vf (or vf/va) is independent of the 
disorder strength vr- However, this is artificial. In the 
present problem, the flow equation of va/vf (or vf/va) 
is derived from the more fundamental equations of va 
and vp, and therefore, is reliable only when va and vf 
both have well-defined fixed points [Ell, H2] ■ If the RG 
equations of va and vp have unphysical values, the run- 
ning equation of Va/vf (or vf/va) becomes meaningless. 
Based on numerical calculations in the presence of ran- 
dom chemical potential, the flow of vf and t>A in vf — va 
space at two representative initial values is presented in 
Fig. ([5]). This implies that the unphysical values of vf 
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FIG. 8: (Color online) Flows of vf,va,va/vf, and vr for 
the t x component in the presence of random gauge poten- 
tial at two representative initial values vao/vfo = 0.1 and 
g = 10~ 3 (The conclusions are independent of these concrete 
values) for: (a) M = t v ; (b) M = r z . Their r z -component 
counterparts are very similar to this case, and thus are not 
presented. 



and va are generated with increasing I. Furthermore, 
they manifest rapid oscillations between positive and un- 
physical negative values as I grows. Therefore, vf and 
va do not reach any stable values due to the interaction 
between nodal QPs and the random chemical potential. 

In order to understand this concretely , we would also 
address a briefly qualitative analysis. The flow equations 
of vp and va in the presence of the random chemical 
potential are 



dvp 
~df 
dvA 



C\ — 6*2 



C\ — Cs 



vh 



2ttvfva J 

v r9 \ 
2itvfva J 



vf, 



VA, 



(76) 
(77) 



where g and vr correspond to the random chemical case. 
In the spirit of RG analysis [l9| , we can obtain the possi- 
ble fixed pints of fermion velocities va and vf by requir- 
ing that 



dvF 
~df 
dvA 
~df 



(Ci - C 2 )v F 
(Cx - C 3 )v A 



2ttva 

vh 

2itvf 



(78) 
(79) 



We assume that v* F and v A correspond to the fixed points. 
If both v* F and v* A are finite, then the above equations 



imply that (Ci — C2)v* A < 
not be satisfied since v% 



7^0. 



(Ci 
If v 



A 



C 3 )v 
= 0, then 



A u Fi 



which can 



vh 



2ttv* a (Cx - C 2 ) 



(80) 



From the expressions for C\ and C 2 , this implies that 
1 cx l/(v A ) 2 , which is clearly inconsistent with the as- 
sumption of v A — 0. Before going to the v* F = case, we 
define Fi — {vf /va)Ci, i = (1,2,3), then the new forms 
of Eqs. (J75J) and ([7H]) become 



(Fx - F 2 )v A 



(Fx-F 3 )^ 



vh 

2ttva ' 

vh 

2itvf 



(81) 
(82) 



If v* F = 0, by both analytical and numerical analysis, we 
found that these equations have no solution. 

In conclusion, the fermion velocities vp and va, as dis- 
cussed above, do not approach any stable values under 
the low energy regime caused by the interaction between 
fermions and random chemical potential as shown in Fig. 
©. Therefore, there is no fixed point of the fermion ve- 
locities vf and va in this case. We interpret this as an 
indicator of the instability of the quantum phase transi- 
tion in the presence of the random chemical potential. 

For completeness, we would like to make brief state- 
ments with various types of disorders present. In general, 
there may be a certain number of kinds of disorders in 
realistic physical problems [22l. l23l. |49|. [53l - [55j . First, we 
consider a case in presence of the random chemical poten- 
tial and random mass (and/or random gauge potential). 
Since the random mass and random gauge potential are 
both irrelevant, the marginal random chemical potential 
dominates and destroys the fixed points (Hf. On the 
other hand, the corresponding fixed points would remain 
for a combination of two irrelevant cases of random mass 
and random gauge potential. Therefore, we can obtain 
the overall effects of disorders by investigating the three 
types of disorders separately. 



VII. VELOCITY RENORMALIZATION 
EFFECTS ON SUPERFLUID DENSITY AND 
CRITICAL TEMPERATURE 

In an actual <i-wave cuprate superconductor, the 
fermion velocities vp and va of the gapless nodal QPs 
are not equal. Indeed, the ratio va/vf may be as low as 
1/20 [1, 5]. As mentioned in Sec. IVI| the value of this 
ratio va/vf will be extremely influenced by A"-ordering 
fluctuations in the proximity of quantum critical points. 
Since the ratio plays an important role by entering a 
number of physical quantities [1], it is natural to expect 
numerous intriguing effects due to the significant changes 
in the ratio va/vf. In this section, we primarily address 
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FIG. 9: (Color online) Flow of (v F ,v A ) for I G [0,600] in the 
presence of random chemical potential at two representative 
initial values «ao/«fo = 0.1 and g — 10~ 3 (The conclusions 
are independent of these concrete values). 



the behavior of supcrfluid density and critical tempera- 
ture at different quantum critical points. 

The superfluid density is an important quantity that 
characterizes the fundamental feature of the d-wave 
superconducting state. In cuprate superconductors, 
the superfluid density is known to exhibit a linear In- 
dependence as @, HH 



p s (T) _ p s (0) 2\n2v F 



7T VA 



T. 



(83) 



with a coefficient proportional to the inverse of the ve- 
locity ratio va/vf- Here, p s (0) — x/a 2 is the super- 
fluid density at zero temperature in the underdoped re- 
gion [l], [56j], where x and a represent the doping con- 
centration and lattice spacing, respectively. For non- 
interacting nodal QPs, the velocity ratio takes a bare 
constant, va as mentioned above. 

An intimately related quantity is the critical tempera- 
ture r c !! m, 



T c oc p s (0) oc — x, 



(84) 



where x is the doping concentration. 

Close to the critical point of the X matrix state, due to 
the strong A-order-parameter fluctuations, the velocity 
ratio flows upon lowering the energy scale. The effects of 
velocity renormalizations should be taken into account. 
To study these in detail, after including the flow of veloc- 
ities, Eqs. (|63l l64l [65]), we can obtain the renormalized 
superfluid density [l3l l38l \El\, 



P%(T)=p s (0)-p n R (T), 



PUT) 



in 



d 2 k 



k B T J (2tt 



(85) 



,(86) 



k B T 



where the velocities uf,a are very complexly dependent 
on k and determined by the running equations in Sec. [V] 
The renormalized critical temperature can be acquired 
by setting 



p s (0) = p R {T c ). 



(87) 



To estimate these more quantitatively, we assume the ul- 
traviolet cutoff A = lOeV, and choose the representative 
bare velocity ratio va/vf ~ 0.1, which is an appropriate 
value for YE^CuaOe+a 

First, we consider the clean limit system. Approaching 
the critical point of the X matrix state, the velocity ra- 
tio flows to va/vf = [20], va/vf — 1, and vp/va = 0, 
corresponding to M = t x , t v , and t z , respectively. This 
leads to the suppression of the superfluid density in case 
the M = t x [38| and growth in the rest of the cases. 
The numerical results for superfluid densities and criti- 
cal temperatures are shown in Fig. (fTU|) . As mentioned, 
based on the fixed point va/vf — > for M = t x , the 
supcrfluid density and critical temperature are both sup- 
pressed. On the other hand, they are slightly increased in 
the case of M = t v . Furthermore, to the case M = t z , it 
exhibits substantial enhancements owing to the running 
flow to vf/va 0. Since plenty of approximations are 
unavoidably employed during the calculations, we need to 
point out here that the quantitative enhancements may 
not be reliable. However, the qualitative increments of 
supcrfluid density and critical temperature would be un- 
ambiguous. 

Next, we discuss the effects caused by disorders, in- 
cluding random mass, random gauge potential, and ran- 
dom chemical potential. According to the analysis in 
Sec. IVI Bl the disorder strength vr is irrelevant both for 
random mass and for random gauge potential. There- 
fore, the behavior of the superfluid density and critical 
temperature mentioned in the previous paragraph are un- 
changed in the presence of these two types of disorders. 
By comparison, quantum phase transitions are unstable 
against the random chemical potential. This suggests 
that the corresponding fixed points do not exist. 

In brief, we conclude that distinct behaviors of physical 
quantities, such as superfluid density and critical tem- 
perature, will be generated by the influence of singular 
velocity renormalization near the presumable A-ordering 
quantum critical points. As a consequence, this leads to a 
possible way to determine the A-ordering quantum criti- 
cal points by seeking these singular behaviors. This may 
be of great help to understand the complicated phase dia- 
gram and many anomalous features of high-temperature 
superconductors. 



VIII. SUMMARY 

In summary, we investigate unusual renormalizations 
of the Fermi velocity vp and gap velocity va of nodal 
QPs close to diverse quantum critical points classified 
by vertex matrices in the d-wave superconducting 
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FIG. 10: (Color online) Superfluid densities [in units of 
zero-temperature superfluid density p s (0)] are calculated for 
A/T c = 10 4 and vao/vfo = 0.1. p a {T) denotes the super- 
fluid density in the absence of effects of quantum fluctua- 
tions. p S ft X (T), Pn V {T), and p^ z (T) correspond to the renor- 
malized superfluid density for M = t x (nematic), M = r y , 
and M — r z , respectively. The critical temperatures labeled 
T x ' y,z were obtained directly from the intersections of the 
curves of superfluid densities with the horizontal axis. It is 
obvious that both the superfluid density and the critical tem- 
perature for M = t x are suppressed. However, their counter- 
parts for M — t v receive slight enhancements, and those for 
M = t z , significant enhancements. 

dome. According to the spirit of the RG method [T9i — 
l2lj . the series of RG equations — both in clean limit 
and in the presence of three types of static disorders, 
namely, random mass, random gauge potential, and ran- 
dom chemical — are derived to one-loop level. A detailed 
RG analysis and numerical computations are given. 

In clean-limit systems, these RG equations lead to 
three distinct fixed points of the velocity ratio: va/vf — > 
20], va/vf — !• 1, and va/vf — > oo for M — t x , t v , and 
r z , respectively. Since the velocity ratio, va/vf, enters 
various physical quantities, these fixed points result in 
abundantly physical observable effects. In particular, we 
calculate the superfluid density and critical temperature. 
They are affected sensitively by approaching these fixed 
points. Superfluid densities and critical temperatures are 
both suppressed, slightly increased, and significantly en- 
hanced for cases M = t x , t v , and r z , respectively. 



Three types of disorder effects on these fixed points are 
also examined. We present a series of coupled RG equa- 
tions of Fermi velocity vf, gap velocity i>a> and disorder 
strength parameter vr . After both analytical and numer- 
ical computations, we find that fixed points obtained in 
the clean limit are robustly stable in the presence of ran- 
dom mass and random gauge potential. However, these 
fixed points are destroyed in the presence of a random 
chemical potential which is marginal. Therefore, this is 
responsible for the instabilities of quantum phase tran- 
sitions. The corresponding effects on physical quantities 
are also studied. Compared to the clean limit, the be- 
haviors of physical quantities controlled by fixed points, 
such as the superfluid density and critical temperature 
depicted in Fig. (|10p . are unchanged when a random 
mass or random gauge potential is included. On the 
other hand, the effects of quantum critical fluctuations 
on physical observables can be neglected due to the fixed 
points broken in the presence of a random chemical po- 
tential. 

The occurrence of quantum phase transition stems 
from competition among ground state phases (5^ - i60j . 
The physical behavior of a system is influenced signifi- 
cantly within a wide scope of the phase diagram, espe- 
cially near quantum critical points where fluctuations are 
divergent [58j . By studying various quantum phase tran- 
sitions in a superconducting dome, the critical behavior 
of physical quantities can be captured. Remarkably, this 
provides a helpful clue to confirm or even locate the very 
existence of quantum critical points by means of detect- 
ing these singular behaviors. Therefore, it will be in- 
structive to understand the complicated phase diagram 
and a number of anomalous features of high-temperature 
superconductors . 
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